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We show that the ground state of the periodic long range Josephson array frustrated by magnetic 
field is a glass for a sufficiently large Josephson energies despite the absence of a quenched disorder. 
Like superconductors, this glass state has non-zero phase stiffness and Meissner response; for smaller 
1-^ , Josephson energies the glass "melts" and the ground state loses the phase stiffness and becomes 

insulating. We find the critical scaling behavior near this quantum phase transition: the excitation 
gap vanishes as (J — Jc)^, the frequency-dependent magnetic susceptibility behaves as xi'-^) oc 

pH . I. INTRODUCTION 

c/3 ' Glass pfra-mation in the absence of intrinsic disorder is a long-standing problem but the last years saw a rapid 

Zj , progresstjlj in the qualitative understanding of this phenomena. Mostly this progress is due to the solutions of 
periodic models that assume a mapping between the periodic model and the appropriate random modelou. The 
j^ ' validity of this assumption is still an open question in a general case but it was shown that at least one periodic model 
allows a direct study of a phase transitiorfl and non-ergodic behavior below the transitioiiO without any reference 
to a disordered model. This model describes a long-range Josephson array in_a^magnetic field and another reason 
''O ' for the interest in this model is that it that can be realized experimentally (cf.Ha for the discussion of experimental 
G ' conditions). 

^ All these results were obtained in the framework of classical statistical mechanics, the glass formation in a regular 

quantum systems has not been addressed. The goal of this paper is ta,fiIL this gap. The problem of a glass formation 
in disordered quantum systems was discussed in a number of workgoliiJ; these works |-s±udied the critical behavior 
near the quantum vitrification transitioncllllj and the properties of the glassy phase itselO using the replica approach. 
They found that the glass phase transition at T = indeed exists; further, it strongly resembles classical (high T) 
phase transition in the same system: the main difference is in the critical exponent of correlation function which 
decays faster than at the classical critical point: D{t) = {Sj{0)Sj{t)) ~ t~^ at T = (cf. D{t) ~ t~-^/^ at non-zero 
p^l T). A surprising result stated irO is that at zero temperature no replica-symmetry-breaking (RSB) is needed for the 

f^ description of the glassy state, i.e. replica-symmetrical solution is stable at T = 0. Since, usually, RSB is believed to 
On , be a signature of non-ergodicity, this result means either absence of non-ergodic behavior at T = or violation of the 
^N . usual relation between RSB and non-ergodicity. We feel that in order to clarify this important question, an approach 
that is free from the ambiguities of the replica method should be employed. 

Understanding of a quantum glass formation in a system with regular Hamiltonian is important for the general 
problem of Quantum Computatioi£3. The reason is that quantum computer is also a quantum system with exponential 
number of states and the process of computation can be viewed as an almost adiabatic change of the external 
fH parameters resulting in a different state. The crucial question is what are the conditions so that such process does not 

Q J lead to the collapse of the density matrix due to the coupling to the environment. This question can be addressed to 
O ■ the spin glass system as well and one can learn about decoherence in a generic large system with exponential number 
^ ' of states from the answer to it. 

-Here wc study the quantum version the long-range Josephson array in a frustrating magnetic field that was suggested 
ina^lll. We consider here only the problem of glass formation, approaching the glass from the "liquid" (i.e. insulating) 
3 ■ side. We show that the quantum version of this problem is described by the same dynamic equations as the quantum 
■ - ' disordered p-spin model studied ir£3. Thus, we explicitly prove that this frustrated quantum system can be mapped 
onto the quantum disordered system in a complete analogy with the situation for classical problems. Further, we 
provide a direct numerical proof that the transition in this model is indeed continuous as conjectured irO and we 
calculate the anomaly of the diamagnetic response associated with this transition. 

Anotlie]p-|and more physicalp-justification of the model is the following. It is well established, both experimentally 
(cf. e.gJl^'Ej) and theoreticallyta that usual nearest-neighbors Josephson arrays made of small superconductive islands 
demonstrate zero-T superconductor-insulator transition as the ratio of the Josephson coupling Ej between the super- 
conductive islands to the Coulomb energy cost Ec = (2e)^/2C for the transfer of the Cooper pair between the islands 
decreases. At small values oi x = Ej/Ec the ground state is an insulator with nonzero Coulomb gap in the excitation 
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spectrum. At nearly-critical values of a; « Xcr the transition between insulating and superconductive states can be 
triggered by application rof a weak magnetic field, producing frustration of the Josepson interaction; moreover, this 
transition can be splitedEj into the sequence of two different transitions: superconductor-metal-insulator. Although 
main qualitative features of these phenomena are understood, there is still no quantitative theory which describes 
quantum phase transitions in 2-dimensional short-range systems, especially in the presence of frustration. Therefore, 
in our attempt to study the origin of a quantum glass state, we have to turn to the simplest (theoretically) model of a 
Josephson array with long-range interaction, which consists of long superconducting wires (instead of islands). It will 
allow us to employ some version of mean-field-theory and reduce the problem to a zero-dimensional quantum theory 
with the interaction that is non-local in time. 

The system of our study is a stack of two mutually perpendicular sets of N parallel thin superconducting wires with 
Josephson junctions at each node that is placed in an external transverse field H . Macroscopic quantum variables 
of this array are the 2N superconducting phases associated with each wire (e.g. the value of the phase of the 
superconducting order parameter at the center of each wire); we will always assume that excitations within individual 
wire can be neglected, so the whole wire is characterized by one phase, (j)m- In the absence of an external field the 
phase differences would be zero at each junction, but this is not possible for finite H , so the phases are frustrated. Here 
we assume that the Josephson currents are sufficiently small so that the induced fields are nejgligible in comparison 
with H (this imposes an important constraint for the experimental realization of this networkQ) . Therefore the array 
is described by the Hamiltonian 

H^nj+nc = -EjyCOs{<Pr.-<Pm-'^ I Ad\) + i^ V C ,;/„ -f- -^ (1) 



^ ^^" ^" he J ' 2 ^^™'", 

where Hj and He represent, correspondingly, Josephson and Coulomb parts of the Hamiltonian, and C„ 
matrix of the capacitances. There are several different contributions to C: self-capacitances of the wires C/ (with 
respect to substrate), the contact capacitances Cj and mutual capacitances of wires C'li- Below we will assume, 
that the self-capacitance is the largest of all, C/ ^ Cu , NCj (the factor N accounts for the fact that there are N 
contacts along each wire). These conditions allow as to neglect all mutual capacitances and consider the matrix 
Cm,n to be diagonal with eigenvalues C/. It is convenient to rewrite the Hamiltonian in terms of "spin" variables 
Sm = e*"^'" . Choosing the Landau gauge for the vector potential and introducing Jq by Ej = -j= so that the transition 
temperature remains constant in the limit A^ — > oo at fixed Jq we get 

^^ Ec 

H = — / ^ S^JmnSn H ^ / , Qn l^j 

ni , n n 

where Qn = —id/d(j)n is the charge operator conjugated to the phase (j>n, Eq — ^, and J^mn is the coupling matrix 

with Jjk = -j= exp{2niajk/N) and 1 < (j, k)<N where j{k) is the index of the horizontal (vertical) wires; Sm = e*"*™ 

where the (pm are the superconducting phases of the 2N wires. Here a = NHP/^q is the flux per unit strip, I is the 
inter-node spacing and $o is the flux quantum. 

Because every horizontal (vertical) wire is linked to every vertical (horizontal) wire, the conujectivity in this model 
is high ( A^) and it is accessible to a mean- field treatment (its classical version was developed intSlj) . For -^ <$^ a < 1 
there exists an extensive number of metastable solutions which minimize the Josephson ("potential") part of the 
Hamiltonian (g); these minima are separated by the barriers that scalea witli A^. A similar (classical) long-range 
network with disorder was previously found to display a spin glass transitionllj for a 3> ;^ ; in the ^sence of short- 
range phase coherence between wires (a ^ 1) it is equivalent to the Sherrington-Kirkpatrick modekEj Physically this 
glassy behavior occurs because the phase differences associated with the couplings, Jjk, acquire random values and 
fill the interval (0, 27r) uniformly. For the periodic case, this condition is satisfied in the "incommensurate window" 
-^ <C a < Ipfor which the magnetic unit cell is larger than the system size so that the simple "crystalline" phase is 
inaccessible. □ There are thus no special field values for which the number of minima of the potential energy are not 
extensive, in contrast to the situation for a > 1. Below we will consider the case 1/A^ <C a <C 1 only. As follows from 
the previous studiesQLl, the characteristic energy scale related to the potential energy Tij is of the order of the glass 
transition temperature of the classical system, Tg ~ Jo/y/a. The zero-T transition we study here is driven by the 
competition between Josephson and Coulomb energies, the scale of the latter being Ec — 4e^/C;. Thus, we expect 



that the quantum transition occurs at Jq/ ^/a ~ Ec- Our goal is to show that such a (continuous) phase transition 
indeed occurs and to study the critical behavior near the transition point. Below, in the main part of the paper, 
we measure all energies in units of £'c, and return to the physical units only in the final expression for the critical 
behavior of the ac diamagnetic susceptibility. 

II. QUANTUM LOCATOR EXPANSION 

We are sping to develop a diagram technique for the Hamitlonian ^ which will be very similar to the one employed 
previouslyu for the classical Langevin dynamics of the same array. The idea is to treat Coulomb part of the Hamiltonian 
as the zero-level approximation, and construct expansion in powers of Josephson coupling constant Jq, keeping all the 
terms of the lowest order in the coctrdination number 1/A^. Thus our approach can be considered as quantum version 
of the Thouless-Anderson-Palmerll3 method. 

The diagram technique for the Matsubara Green function 

G^A^) - -(T,s™(t)s],(0)), s{t) = e--^se^^ (4) 

is closely related to the one developed inEl. Dyson equation for the frequency-dependent matrix Green function reads 
(note that in our units Ec = 1): 

G^ = -^- ^ ^ (5) 

Gji - (JJt)G^ 

where we introduced the local Green functions G^ that is irreducible with respect to the Jy lines. The matrix {J J^)ij 
depends only on the "distance" i — j and acquires a simple form in Fourier space {J ,P)p = {Jq /a)6{a7T— \p\); therefore 
in this representation 



(aTT— IpI) 9{\p\ — an) 
Gji - 4g^ G-^ 



gup)- :' ,r + ^^'^L~""^ (6) 



Diagrammatically Eq.(|g) and the equation for the irreducible function G^j are represented by the graphs shown in 
Fig. 1. 
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Note that the equation for G is written in the lowest nontrivial order in a. Indeed, it is seen from Eq.(|6|) that nontrivial 
part of the Green function which contains critical slowing down, is of relatively small weight ~ a. It is this long-time 
part of G(j which enters 3-line diagram on Fig.l and makes it proportional to a^] more complicated diagrams either 
contain even higher powers of a, or are small as 1/N . Since the second diagram on Fig.l contains single-site functions 
only, the whole system of equations can be written in the form 

G{u:) = {l-a)G{u) + G{u); G{u^) ^ l^t\ / (^) 

1 - jQG^[uj)/a 

G{uj) = Gq{uj) + Y.{oj)- E(cj) ^ ( ^ j xl j G^{t) exp(iwt) dt (8) 

Here X3 ~ Ij as inti, is a static value of four-point vertex denoted as a square box in Fig.l (we assume that, like inH, 
the main critical anomaly is contained in the 2-point Green function alone). Equations (|3,§) should be solved with 
obvious initial condition: 

G{t^O)^ f^G{cj)^l (9) 



Similar normalization condition in the classical problem was sufBcient to determine G{liJ = 0) exactlyul. The same 
calculation seems to be difficult for the present quantum problem and we will not perform it here. Instead, we will use 
general properties of the function 00(0;), namely: i) G(j{Q) ^ 1, and ii) Go{uj) is analitic at low uj and has characteristic 
frequency scale of the order of 1. In doing so, we will not determine exact position of the phase transition (i.e. critical 
value Jc of the coupling strength Jq), but we will show the existence of continuous transition and find the form of 
critical scaling. 

Let us first analyze equations (0,||,||) omitting the term with S, and using simplest interpolation Go{i-^) = {X+uj"^)'^. 
Then initial condition (9|) gives us the equation for A: 

- ' +T^^ (10) 
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where g — Jo/y/a. Thus A ~ 1 as long as g < 1. On the other hand, at 17 ^ 1 the solution is A — g = a ~ (a/4)^. 
The value of a determines the asymptotic decay rate of the Green function 

G(i) = ^exp(-|t|V^) (11) 

with S being neglected. It will be seen below that a '^ a and thus A ~ 1 near the phase transition point g = gc (we 
will see also that S ~ a and thus it is much smaller than the oj^ term at high frequencies uj 3> a^"). It means that 
the parameter a can be considered as a smooth function of g in the vicinity of gc- Clearly, this conclusion does not 
depend on the model form of Gq{uj) used in the above analysis. 

Now we re-introduce S(w) into the equations for G{uj) and focus on its low-frequency behavior at w < -./a: 

G{i^) = ^^F^r-i 7 , (12) 

Y,{uj) ^g^ f G^{t) exp{iu;t) dt (13) 

where g = gx^ ^ g- Strictly speaking, Eqs.( |l2[|l3| ) do not form closed system since a should be determined 
with the use of Eq.dS) which contains high-frequency contributions. However, in this high-frequency region (which 
produces the main contribution to the normalization condition (^ ) the contribution of S(tt') can be neglected and 
thus a can be treated as an external control parameter which governs the transition. Green function defined by the 
Eqs.( [12| , [l3| ) acquires singularity when 2I](0) — a. To find the form of this singularity, we make use of the scaling 
Anzats G{t) = qt~^ and neglect w^ term in the denumenator of Eq.(^2|). Then we find v = 1/2 and q ^ g~^o^l'^ . This 
critical-point solution matches the short-time asymptotics (|ll|) at t ~ a~^l'^ . The estimate for E(0) which follows 
from the above scaling Anzats, 

00 
I](a. = 0)^.gV/^-5V«^/^ 

gives S(0) « a at g ^ 1 and a ~ a, as it was expected. These estimates show that second-order phase transition 
with critical slowing down may indeed occur in the above range of parameters. In the next section we will study the 
vicinity of the critical point in more details. 

III. GREEN FUNCTION NEAR THE T = TRANSITION POINT 

To study the form of the critical singularity, it is convinient to define universal scaling functions Q[ui) and <j{lo) 
which do not contain small parameter q <C 1, and the parameter h measuring the proximity to the critical point: 

G{u:)^Q(i::Y aS(w) = cr(a)); (i = w/^a; 6 = (a - 2E(0))/a (14) 

Now Eqs.( |l^ , [l^ ) acquires the following form: 

^(^) = ,^2(a(0)-^(^)) ^ "^^^ = ^' / ^'^^~^ "^^"^ '' ^''^ 



Exactly at the critical point 6 = the solution of Eq. ( |15| ) is 

'7r\i/4 



Qi^) = ( g 



r^/2p|-i/2 



(16) 



Consider now the vicinity of the critical point, < fe ^ 1. It is clear form the form of the solution (|16| ) that the 
similar result should be valid at a) 3> 6^. Next we focus on the long-time, low-w region, cj <C 6^ and will look for the 
purely exponential solution 

g{i) = giexpi-i/To). (17) 

This type of asymptotic solution is known to exist in the classical version of the same model (cf-O'Ll) . In the present 
problem, one can show, considering analytic structure of Eqs.(p^), that at 6 > the singularity of Q{u)) which is 
closest to the real axis of cj, is necessarily a simple pole at some w = i/TQ] the next singularity may exist at u) > Si/tq. 
Solution of Eqs.(Q3) with the Anzats (O) in the region i^ tq determines parameters tq and Gi as functions of b: 

(32/27)1/2 1 



To 
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(27/2)i/2g3 5 



(18) 



This solution is similar to the one found inti; however, an important difference is that in the present case the prefactor 
Gi scales to zero at the critical point 5 = 0. 



The full solution in the vicinity of the transition point should contain both (16) and ( |18| ) as asymptotic solutions, 
and can be written in the form 
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Gi exp 



t 
To 



(19) 



where f{x) is some scaling function approaching constant at a; = and fast decaying at x — > cx), and ri < tq/S. To 
confirm an existence of this type of solution, we solved Eqs.([l5|) numerically for several values of 6 ^ 1. The results of 
this computation are shown on Fig. 1. Clearly, all three functions Gii^) coincide in the high-o; region, there they are 
close the square-root asymptotic dlq). Low-frequency parts (for ui < 0.08) of these solutions can be made coinciding 
by a proper rescaling of their arguments, to* = Auj. Fig. 2 demonstrates linear relation between b^^ and the scaling 
coefficient A, as it was suggested by Eqs.( |l^Jl^ ). 
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FIG. 1. Low-frequency asymptotic behavior of G{u>) at different 6 at T = 0. 
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FIG. 2. The relation between scaling parameter A and the proximity to the transition point 6 at T = 0. 
These results confirm the existence of the T = critical behavior of the type of Eq. dl^) . 



IV. CRITICAL BEHAVIOR AT T > 

The above results refer to the zero-T phase transition controlled by the single parameter g — Jo/^/a. We found 
that this phase transition is a continuouSj-pne and the corresponding critical behavior differs considerably from the 
one found in an analogous classical modelQ. In particular, at T = critical point g — gc there is no "plato" solution 
with approximately constant G{t) at t — > oo which is known to be a peculiar property of regular classical glasses. Now 
we consider low but non-zero temperatures T = f3~^ and will find how "classical" critical scaling "grow up" from the 
"quantum" background; we also find the low-temperature shape of the phase transition line on the plane {T,g). 

Green function is defined now at discrete frequencies a;„ — 2TrnT and the equations (|l2|), ( |l3[ ) can be written as 



GK) = 



a — 2Y,{uJn) + w^ 



EM(t^n) = g / G^(t)exp(iw„i) dt 





(20) 



It will be convenient now to perform analytic continuation of Eqs.(pOD and rewrite them in terms of real-time correlation 
function D{t) = {[S{t), S{0)] + ) and response function x(t) = i[S{t),S{0)]-e{t). The functions G(w„), D{uj) and xi^) 
are related as follows: 



G{-iu; + 7]) = x(t^), V^ +0 , D{uj) = Imx(w) coth(cj/2r) 

After analytic continuation the Eqs.(|2C|) can be written as 

x(^) = ~ 



5-2E(t^)' 



S(w) = Sg" / D\t)x{t) (cxpiiLut) - 1) dt; a^a- 2I](w = 0) 



(21) 



(22) 



where we omitted ui^ term which is irrelevant in the vicinity of the critical point. Equations (|2|) form (together with 
the Fluctuation-Dissipation relation (second of Eqs. (|2l|)) a closed set which determines critical singularity at T > 0. 
Formally Eqs.(p^ coincide with the corresponding "classical" equations fromO, the only difference is in the form of 
the Fluctuation-Dissipation relation. 

Let us consider low temperature region T <^ ^/a. As long as we are interested in the long-time asymptotic 
t ^ 1/T, the correlation and response functions are related by classical FDT: D{u!) ~ 2T/uj Imx(a;). Characteristic 



times, which is relevant in (22), are also belong to classical region t ^ 1/T. Therefore the correlation function at 
the transition point has the same critical behavior as in the classical case: limt^oo D{t) = q. However parameter 
a = A — 5 is determined by the "quantum" region of frequencies uj ^ T, i.e. by the equation (Rw. The substitution 



of this expression to the equations (E3) allows us to find 



q ^ a^l^T^'^- ~a ~ a3/4j.i/2 (23) 

In the short-time domain t <C T~^ the zero-T critical solution with D{t) ^ Q^l'^t~^l'^ is valid. Equation (p 
demonstrates the way the "classical" solution with nonzero linii^oo ^(0 grows up with the temperature increase 

V. DIAMAGNETIC RESPONSE NEAR THE TRANSITION POINT 

Correlation and response functions -D(i) and x(i) are not directly measurable in our system, but they can be used 
in order to calculate measurable physicaLquantity which is dynamic diamagnetic susceptibility xa<('^)j like it was 
done previously for the classical problemEl. Total magnetic moment induced by time-dependent external magnetic 
field is given by 

^ ^ mn 

where Srnn — iTnnJmn^- Then magnetic susceptibihty xm c^n. be found making use of the Kubo formula: x_a/i(^— ^') = 
i [M{t), M^t')] 6{t) which leads to the expression 

Xm{uj) = f ^ j I' J (e'"* - 1) ReTrJx{t)JD{t)dt . (25) 

Here we omit the term, containing irreducible four-spin correlator (of the order of 1/N), and take inta.account that 
M.{H = 0) = 0. Note, that equation (|25| ) formally coincides with classical formula for magnetic responses. The matrix 
functions D{t) and xit) contain elements (denoted by superscript *^°^) belonging to the same (horizontal or vertical) 
sublattice of our array, as well as "offdiagonal" elements (with superscript '"'^^) which describe correlation of phases on 
wires of different type (horizontal/ vertical). Relation between these functions is as follows: x(w) — JG{uj)x{ujy ' . 
Thus, the expression for magnetic susceptibility has the following form: 



XA.(-) = (^)7^V^'4/M (26) 



where 

I{u) = I \5{t- h) - '""^" ^'' 1 x{ti)D{h) (e'"* - 1) e{t)dtdh (27) 



J^G^t — ti) 



Near the transition point only long-time parts of all the function in ( p7| ) are relevant, and this expression can be 
reduced to the form 

/(u;) = (S(c.) - E(0)) / x{t)D{t)e'-'dt , (28) 



where the first factor came from the first brackets in (M); note that it vanishes in the limit uj 
Using the solution (^^ we obtain at the quantum critical point J = Jc- 

X 1/4 



^M = ^(^) V^^lii'^ (29) 



Near the T = transition point at high enough frequencies uj ^ {J/Jc ~ 1)^Q; ■^/^ equation ( p9| ) still holds. In 
opposite case of low frequencies 

'("^ = 81(J. - J)3 g ■ (30) 

Note that the parameter g (which is known up to the factors of order 1 only) does not enter the low-w asymptotic of 

Making use of the Eqs.( p6| , p9| , p0[ ) and returning to the original units of frequency, we obtain finally ac diamagnetic 
susceptibility near the quantum transition point: 

XmH - (^^'l'N^^-:^^^V^^JcJe^HuCile\ c » ^'^^^,~Jf (31) 

2e \ V ^' \ \.5 2Qa^/2 J^ ^ ^^ Ci (J - J,) 



The above expressions are valid at the frequencies oj ^ T/Ti^ otherwise "islassical" asymptotic for the Green functions 
should be used and will lead to the frequency dependencies like those inB. 

VI. CONCLUSIONS 

We have shown that regularly frustrated long-range Josepshon array has a quantum (zero-temperature) phase 
transition between Coulomb-dominated insulator phase and a superconductive state. This transition happens when 
the nearest-neighbors Josephson coupling exceeds the critical value, Jy ^ N~'^/'^^/ae^ /Ci, where C/ is the self- 
capacitance of an individual wire. 

We found that quantum critical behavior of the model at J — > Jc is different from that of an analogous classical 
systemEl: at the quantum critical point D{t) ~ t~^/'^ while at the classical critical point q = limj^oo D{t). However, 
at any non-zero temperature a "classical" type of asymptotic behavior is recovered at the longest-times, t 3> h/T, 
leading to g ex T^/^. Near the T = critical point the gap in the excitation spectrum decreases as t^^ cx {Jc — J)^- 

Near the phase transition the effective inductance £ of the array (defined by £ c>c (£^m{'^) / dLiP'\ui-t{)) diverges as 
{Jc — J)^^ ^ this shows that the glass state has a macroscopic phase rigidity (cf. alsoEj). Right at the critical point 
we find unusual frequency behavior of the complex diamagnetic susceptibility, xm{^) oc A/iujhiu}. 

Frustrated nature of couplings in our array and comparison with the previous resultsQ on the classical version of 
the same model indicates that the high- J state is a quantum glassy superconductor. The T = nonergodic properties 
(irreversibility, ageing) remain an open question; note here that recent studyllj of nonequilibrium glassy behavior 
in a p-spin spherical quantum model assumed strongly dissipative (overdamped) dynamics of the whereas dynamics 
relevant for the Josephson array at T = must be underdamped. 
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